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On the Fundamental Formulae of Dynamics. 

By J. W. Gtibbs, New Haven, Conn. 

Formation of a new Indeterminate Formula of Motion by the Substitution of the 

Variations of the Components of Acceleration for the Variations 

of the Coordinates in the usual Formula. 

The laws of motion are frequently expressed by an equation of the form 

(1) X l(X- mx) Sx + ( Y— my) hy + {Z— mz) 8z] = , 

in which 

n\ denotes the mass of a particle of the system considered, 
x, y, z its rectangular coordinates, 
x, y, z the second differential coefficients of the coordinates with respect to 

the time, 
X, Y, Z the components of the forces acting on the particle, 
Sx, 8y, Sz any arbitrary variations of the coordinates which are simultane- 
ously possible, and 
2 a summation with respect to all the particles of the system. 
It is evident that we may substitute for Sx, hy, hz any other expressions 
which are capable of the same and only of the same sets of simultaneous 
values. 

Now if the nature of the system is such that certain functions A, B, etc. 
of the coordinates must be constant, or given functions of the time, we have 
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These are the equations of condition, to which the variations in the general 
equation of motion (1) are subject. But if A is constant or a determined 
function of the time, the same must be true of A and A. Now 

/ _, /dA • dA ■ , dA -\ 
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and 

v v fd A ■■ dA- dA-\ rj. 

where H represents terms containing only the second differential coefficients 
of A with respect to the coordinates, and the first differential coefficients of 
the coordinates with respect to the time. Therefore, if we conceive of a varia- 
tion affecting the accelerations of the particles at the time considered, but not 
their positions or velocities, we have 

(dA .•• , dA ,•• , dA 



(3) 



*2 = s(£tf + £* + £S) = o, 



and, in like manner, 



«= 2 (S^+f^+S^)=°. 



etc. 

Comparing these equations with (2), we see that when the accelerations 
of the particles are regarded' as subject to the variation denoted by h, but not 
their positions or velocities, the possible values of hx, hy, hz are subject to pre- 
cisely the same restrictions as the values of hx, hy, hz, when the positions of 
the particles are regarded as variable. We may, therefore, write for the 
general equation of motion 

(4) 2 l(X - mx) hx + ( Y- my) hy + (Z— mz) hz] = , 

regarding the positions and velocities of the particles as unaffected by the 

variation denoted by 8, — a condition which may be expressed by the equations 

8x = 0, $y = 0, $z = 0, 
Sx = 0, Sy = 0, 8z = 0. 

We have so far supposed that the conditions which restrict the possible 
motions of the systems may be expressed by equations between the coordinates 
alone or the coordinates and the time. To extend the formula of motion to 
cases in which the conditions are expressed by the characters < or > , we 
may write 
(6) s\{X—mx) hx + (Y— my) h'y + {Z - mz) &s] =0. 

The conditions which determine the possible values of Sx, hy, hz will not, 
in such cases, be entirely similar to those which determine the possible values 
of hx, hy, hz, when the coordinates are regarded as variable. Nevertheless, 



(5){ 
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the laws of motion are correctly expressed by the formula (6), while the 
form ill a 

(7) X Ux-mx) hx + (Y-my) hy + (Z-m'z) hz~\ = 0, 

does not, as naturally interpreted, give so complete and accurate an expres- 
sion of the laws of motion. 

This may be illustrated by a simple example. 

Let it be required to find the acceleration of a material point, which, at 
a given instant, is moving with given velocity on the frictionless surface of a 
body (which it cannot penetrate, but which it may leave), and is acted on by 
given forces.- For simplicity, we may suppose that the normal to the surface, 
drawn outward from the moving point at the moment considered, is pai'allel 
to the axis of X and in the positive direction. The only restriction on the 
values of hx, hy, hz is that 



Formula (7) will therefore give 



hx = 0. 



■■> X ■■ Y ■■ Z 

x — — , y = — ■, z— — ■. 

The condition that the point shall not penetrate the body gives another 
condition for the value of x. If the point remains upon the surface, x must 
have a certain value N, determined by the form of the surface and the velocity 
of the point. If the value of x is less that this, the point must penetrate the 
body. Therefore, 

x^ If. 
But this does not suffice to determine the acceleration of the point. 

Let us now apply formula (6) to the same problem. Since x cannot be 
less than N, 

\f'x = N, hx = 0. 
This is the only restriction on the value of hx, for if x > If, the value of hx is 
entirely arbitrary. Formula (6), therefore, requires that 

> X 



if x =z N, x 



m 



X 
but if x > N, x — — : 

m 

— that is, (since x cannot be less than If) , that x shall be equal to the greater 
of the quantities iVand — , or to both, if they are equal, — and that 
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The values of x, y, z are therefore entirely determined by this formula in 
connection with the conditions afforded by the constraints of the system.* 

The following considerations will show that what is true in this case is 
also true in general, when the conditions to which the system is subject are 
such that certain functions of the coordinates cannot exceed certain limits, 
either constant or variable with the time. If certain values of hx, hy, hz (with 
unvaried values of x, y, z, and x, y, z) are simultaneously possible at a given 
instant, equal or proportional values with the same signs, must be possible for 
hx, hy, hz immediately after the instant considered, and must satisfy formula 
(1), and therefore (6), in connection with the values of x, y, z, X, T, Z imme- 
diately after that instant. The values of x, y, z, thus determined, are of 
course the vei*y quantities which we wish to obtain, since the acceleration of 
a point at a given instant does not denote anything different from its accelera- 
tion immediately after that instant. 

For an example of a somewhat different class of cases, we may suppose 
that in a system, otherwise free, x cannot have a negative value. Such a con- 
dition does not seem to affect the possible values of hx, as naturally inter- 
preted in a dynamical problem. Yet, if we should regard the value of hx in 
(7) as arbitrary, we should obtain 

"— Z 
m 

which might be erroneous. But if we regard hx as expressing a velocity of 
which the system, if at rest, would be capable, (which is not a natural signifi- 
cation of the expression,) we should have hx = 0, which, with (7), gives 

••> X 

m 
This is not incorrect, but it leaves the acceleration undetermined. If we 
should regard hx as denoting such a variation of the velocity as is possible 
for the system when it has its given velocity (this also is not a natural 

*The failure of the formula (7) in this case is rather apparent than real ; for, although the formula appa- 
rently allows to x, at the instant considered, a value exceeding both N and — , it does not allow this for any 
interval, however short. For if x < N, the point will immediately leave the surface, and then the formula 
requires that x = — . 
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signification of the expression), formula (7) would give the correct value 
of x except when x — 0. In this case (which cannot be regarded as excep- 
tional in a problem of this kind), we should have hx = 0, which will leave x 
tmcletermined, as before. 

The application of formula (6), in problems of this kind, presents no 
difficulty. From the condition 

x = 0, 

we obtain, first, if x = , x = , 

then, if x = and x = 0, $x = 0, 

which is the only limitation on the value of hx. With this condition, we 

deduce from (6) that either 

x =. 0, x = 0, and x = — ; 

X 

or x = — . 

m 

jf •• X 

That is, if # = 0, x has the greater of the values — and 0; otherwise, x = — . 

° mm 

In cases of this kind also, in which the function which cannot exceed a 
certain value involves the velocities (with or without the coordinates), one 
may easily convince himself that formula (6) is always valid, and always 
sufficient to determine the accelerations with the aid of the conditions afforded 
by the constraints of the system. 

But instead of examining such cases in detail, we shall proceed to con- 
sider the subject from a more general point of view. 



Comparison of the New Formula with tlie Statical Principle of Virtual Velocities.— 
Case of Discontinuous Changes of Velocity. 

Formula (1) has so far served as a point of departure. The general 
validity of this, the received form of the indeterminate equation of motion, 
being assumed, it has been shown that formula (6) will be valid and suffi- 
cient, even in cases in which both (1) and (7) fail. We now proceed to show 
that the statical principle of virtual velocities, when its real signification is 
carefully considered, leads directly to formula (6), or to an analogous formula 
for the determination of the discontinuous changes of velocity, when such 
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occur. This will be the case even if we start with the usual analytical expres- 
sion of the principle 

(8) 2 (XSx + Yhy + Zhz) = 0, 

to which, at first sight, formula (6) appears less closely related than (7). For 
the variations of the coordinates in this formula must be regarded as relating 
to differences between the configuration which the system has at a certain 
time, and which it will continue to have in case of equilibrium, and some 
other configuration which the system might be supposed to have at some sub- 
sequent time. These temporal relations are not indicated explicitl}' in the 
notation, and should not be, since the statical problem does not involve the 
time in any quantitative manner. But in a dynamical problem, in which we 
take account of the time, it is hardly natural to use hoc, hy, Sz in the same 
sense. In any problem in which x, y, z are regarded as functions of the time, 
&r, Sy, Sz are naturally understood to relate to differences between the con- 
figuration which the system has at a certain time, and some other configura- 
tion which it might (conceivably) have had at that time instead of that which 
it actually had. 

Now when we suppose a point to have a certain position, specified by 
x, y, z, at a certain time, its position at that time is no longer a subject of 
hypothesis or of question. It is its future positions which form the subject of 
inquiry. Its position in the immediate future is naturally specified by 

x + xdt + -~- xdt 2 + etc., y + ydt + -^ y^ + etc., z-\- zdt -{- — zdt 2 + etc., 

and we may regard the variations of these expressions as corresponding to 
the Sx, 8y, Sz of the statical problem. It is evidently sufficient to take account 
of the first term of these expressions of which the variation is not zero. Now, 
x, y, z, as has already been said, are to be regarded as constant. With respect 
to the terms containing x, y, z, two cases are to be distinguished, according as 
there is, or is not, a finite change of velocity at the instant considered. 

Let us first consider the most important case, in which there is no discon- 
tinuous change of velocity. In this case, x, y, z are not to be regarded as 
variable (by <$), and the variations of the above expressions are represented by 

-=- Sx dt 2 , — hy df, — hz dt 2 , 

a it Z 
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which are, therefore, to be substituted for &r, $y, hz in the general formula of 
equilibrium (8) to adapt it to the conditions of a dynamical problem. By 

this substitution (in which the common factor — dl 2 may of course be omit- 

ted), and the addition of the terms expressing the reaction against accelera- 
tion, we obtain formula (6). 

But if the circumstances are such that there is (or may be) a discontinuity 
in the values of x, y, z at the instant considered, it is necessary to distinguish 
the values of these expressions before and after the abrupt change. For 
this purpose, we may apply x, y, z to the original values, and denote the 
changed values by x + Ax , y + Ay , z + Az. The value of x at a time 
very shortly subsequent to the instant considered, will be expressed by 
x + (x + Ax) dt -f etc., in which we may regard Ax as subject to the varia- 
tion denoted by S. The variation of the expression is therefore 8Ax dt. In- 
stead of — mx, which expresses the reaction against acceleration, we need in 
the present case — Ax to express the reaction against the abrupt change of 
velocity. A reaction against such a change of velocity is, of course, to be 
regarded as infinite in intensity in comparison with reactions due to accelera- 
tion, and ordinary forces (such as cause acceleration) may be neglected in 
comparison. If, however, we conceive of the system as acted on by impulsive 
forces, (i. e. such as have no finite duration, but are capable of producing 
finite changes of velocity, and are measured numerically by the discontinui- 
ties of velocity which they produce in the unit of mass,) these forces should 
be combined with the reactions due to the discontinuities of velocity in the 
general formula which determines these discontinuities. If the impulsive 
forces are specified by X, Y, Z, the formula will be 

(9) [(X — mAx) h Ax + (Y — mAy) SAy + (Z — mAz) 8Az) = . 

The reader will remark the strict analogy between this formula and (6), 

(aOO Utf C?S 

which would perhaps be more clearly exhibited if we should write -j-, -j~, -j- 

for x, y, z in that formula. 

But these formulae may be established in a much more direct manner. 
For the formula (8), although for many purposes the most convenient expres- 
sion of the principle of virtual velocities, is by no means the most convenient 
for our present purpose. As the usual name of the principle implies, it holds 
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true of velocities as well as of displacements, and is perhaps more simple and 
more evident when thus applied.* 

If we wish to apply the principle, thus understood, to a moving system 
so as to determine whether certain changes of velocity specified by Ax, Ay, Az 
are those which the system will really receive at a given instant, the veloci- 
ties to be multiplied into the forces and reactions in the most simple appli- 
cation of the principle are manifestly such as may be imagined to be com- 
pounded with the assumed velocities, and are therefore properly specified by 
&Ax, $Ay, $Az. The formula (9) may therefore be regarded as the most direct 
application of the principle of virtual velocities to discontinuous changes of 
velocity in a moving system. 

In the case of a system in which there are no discontinuous changes of 
velocity, but which is subject to forces tending to produce accelerations, when 
we wish to determine whether certain accelerations, specified by x, y, z, are 
such as the system will really receive, it is evidently necessary to consider 
whether any possible variation of these accelerations is favored more than it 
is opposed by the forces and reactions of the system. The formula (7) ex- 
presses a criterion of this kind in the most simple and direct manner. If we 
regard a force as a tendency to increase a quantity expressed by x, the pro- 
duct of the force by Sx is the natural measure of the extent to which this 
tendency is satisfied by an arbitrary variation of the accelerations. The prin- 
ciple expressed by the formula may not be very accurately designated by the 
words virtual velocities, but it certainly does not differ from the principle of 
virtual velocities (in the stricter sense of the term), more than this differs from 
that of virtual displacements, — a difference so slight that the distinction of the 
names is rarely insisted upon, and that it is often very difficult to tell which 

* Even in Statics, the principle of virtual velocities, as distinguished from that of virtual displacements, has 
a certain advantage in respect of its evidence. The demonstration of the principle in the first section of the 
Mecanique Analytique, if velocities had been considered instead of displacements, would not have been exposed 
to an objection, which has been expressed by M. Bertrand in the following words : " On a objecte, avec raison, 
a cette assertion de Lagrange I'eznmple d'un point pesant en equilibre au sommot le plus eleve d'une courbe; il 
est evident qu'un deplacement infiniment petit le ferait descendre, et, pourtant, ce deplacement ne se produit 
pas." (Mecanique Analytique, troisieme edition, tome 1, page 22, note de M. Bertrand.) The value of z (the 
height of the point above a horizontal plane) can certainly be diminished by a displacement of the point, but 
value of z is not affected by any velocity given to the point. 

The real difficulty in the consideration of displacements is that they are only possible at a time subsequent 
to that in which the system has the configuration to which the question of equilibrium relates. We may make 
the interval of time infinitely short, but it will always be difficult, in the establishing of fundamental princi- 
ples, to treat a conception of this kind (relating to what is possible after an infinitesimal interval of time) with 
the same rigor as the idea of velocities or accelerations, which, in the cases to which (9) and (6) respectively 
relate, wo may regard as communicated immediately to the system. 
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form of the principle is especially intended, even when the principle is enun- 
ciatedor discussed somewhat at length. 

But, although the formulae (7) and (9) differ so little from the ordinary- 
formulae, they not only have a marked advantage in respect of precision and 
accuracy, but also may be more satisfactory to the mind, in that the changes 
considered (to which S. relates), are not so violently opposed to all the possi- 
bilities of the case as are those which are represented by the variations of the 
coordinates.* Moreover, as we shall see, they naturally lead to various 
impoi'tant laws of motion. 

Transformation of the New Formula. 

Let us now consider some of the transformations of which our general 
formula (7) is capable. If we separate the terms containing the masses of 
the particles from those which contain the forces, we have 

(10) 2 (Xtix + Yly + Zhz) -2 [y m & (x 2 +f + z 2 )] = , 
or, if we write u for the acceleration of a particle, 

(11) 2 {XSx + Yhy + Zhz) —52 (-|- m« 2 ) = . 

If, instead of terms of the form Xhx, or in addition to such terms, equa- 
tion (1) had contained terms of the form Php, in which p denotes any quantity 
determined by the configuration of the system, it is evident that these would 
give terms of the form Php in (7), (10) and (11). For the considerations 
which justified the substitution of Sx, hy, hz for hx, hy, hz in the usual formula 

*It may have seemed to some readers ot the Meeanique Anatytique — a work of which the unity of method 
is one of the most striking characteristics, and that to which its universally recognized artistic merit is in great 
measure due — that the treatment of dynamical problems in that work is not entirely analogous to the treatment 
of statical problems. The statical question, whether a system will remain in equilibrium in a given configura- 
tion, is determined by Lagrange by considering all possible motions of the system and inquiring whether there 
is any reason why the system should take any one of them A similar method in dynamics would be based 
upon a comparison of a proposed motion with all other motions of which the system is capable without violating 
its kinematical conditions. Instead of this, Lagrange virtually reduces the dynamical problem to a statical 
one, and considers, not the possible variations of the proposed motion, but the motions which would be possible 
if the system were at rest. This reduction of a given problem to a simpler one, which has already been solved, 
is a method which has its advantages, but it is not the characteristic method of the Meeanique Analytique. That 
which most distinguishes the plan of this treatise from the usual type is the direct application of the general 
principle to each particular case. 

The point is perhaps of small moment, and may be differently regarded by others, but it is mentioned here 
because it was a feeling of this kind (whether justified or not) and the desire to express the formula of motion 
by means of a maximum or minimum condition, in which the conditions under which the maximum or mini- 
mum subsists should be such as the problem naturally affords, (Gauss's principle of least constraint being at the 
time unknown to the present writer, and the conditions under which the minimum subsists in the principle of 
least action being such that that is hardly satisfactory as a fundamental principle,) which led to the formulae 
proposed in this paper. 
16 
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were in no respect dependent upon the fact that x, y, z denote rectangular 
coordinates, but would apply equally to any other quantities which are deter- 
mined by the configuration of the system. 

Hence, if the moments of all the forces of the system are represented 
by the sum J» {Pdp) > 

the general formula of motion may be written 

(12) » (Php) - SX (1 »m 8 ) = . 
If the forces admit of a force-function V, we have 

8V-S2 (i-nuA = 0, 
or W / 

(13) a[f-2 (-|-mt« 2 )] = 0. 

But if the forces are determined in any way whatever by the configura- 
tion and velocities of the system, with or without the time, X, Y, Z and P 
will be unaffected by the variation denoted by 6, and we may write the 
formula of motion in the form 

(14) 62 (Xx + Yy + Zz - y miA = , 

or 

(15) 5 [# {Pi) -S(|m« 2 )] = 0. 

If the forces are determined by the configuration alone, or the configura- 
tion and the time, &X=Q, 6Y=0, 8Z=0, 8P = 0, and the general formula 
may be written _. 

(16) 5 [J. 2 {Xx +Yy + Zz) -2 (-L rotf)] = 0, 
or 

(17) a[£*(i$)-2(ymtf)] = 

The quantity affected by & in any one of the last five formulae has not 
only a maximum value, but absolutely the greatest value consistent with the 
constraints of the system. This maybe shown in reference to (15) by giving 
to p, x, y, z, contained explicitly or implicitly in the expression affected by <5, 
any possible finite increments p', x", y', z', and subtracting the original value 
of the expression from the value thus modified. Now, 

£[P(ii+iO]-2[!m{(*+iT+(y+^ 

= $ {Pp') — 2 \m {xx' + yj/+z z')] — 2 [y m (sf 1 + f + z' 2 )] 
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But since p', if, y', z' are proportional to and of the same sign with possible 
values of hp, hx, by, hz, we have, by the general formula of motion, 

$ (Pf) — 2 [i» (xx 1 + ytf + zz'j] = . 

The second member of the preceding equation is therefore negative. The 
first member is therefore negative, which proves the proposition with respect 
to (15). The demonstration is precisely the same with respect to (13) and 
(14), which may be regarded as particular cases of (15). 

To show the same with regard to (16) and (17), we have only to observe 
that the quantities affected by <5 in these formulae differ from those affected 
by the same symbol in (14) and (15) only by the terms 

2 (Xx + Yy + Zz) and £ (Pp) , 
which Avill not be affected by any change in the accelerations of the system. 

When the forces are determined by the configuration (with or without 
the time), the principle may be enunciated as follows: The accelerations in 
the system are always such that the acceleration of the rate of work done by 
the forces diminished by one-half the sum of the products of the masses of the 
particles by the squares of their accelerations has the greatest possible value. 

The formula (17), although in appearance less simple than (15), not only 
is more easily enunciated in words, but has the advantage that the quantity 

-j- $ (Pp) is entirely determined by the system with its forces and motions, 

which is not the case with $ (Pp)- The value of the latter expression depends 
upon the manner in which we choose to represent the forces. For example, 
if a materia] point is revolving in a circle under the influence of a central 
force, we may write either Xx -f- Yy + Zz or Br for Pp, B and r denoting 
respectively the force and radius vector. Now Xx + Yy + Zz is manifestly 

unequal to Br. But Xx + Yy + Zz is equal to Br, and — (Xx + Yy + Zz) 

at 

is equal to -=- (Br). 

It may not be without interest to see what shape our general formulae 
will take in one of the most important cases of forces dependent upon the 
velocities. If a body which can be treated as a point is moving in a medium 
which presents a resistance expressed by any function of the velocity, the 
terms due to that resistance in the general formula of motion may be expressed 
in the foi-m 

* [V(») ~' x + <?> («) ~-'y + * 00 7 *] > 
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where v denotes the velocity and <?> (v) the resistance. But 

xx , yy , zz dv 

v v v at 

The terms due to the resistance reduce, therefore, to 

3 O («) v] , 

where /' denotes the primitive of the function denoted by <?>. 

Discontinuous Changes of Velocity. — Formula (9), which relates to discon- 
tinuous changes of velocity, is capable of similar transformations. If we set 

w 1 = Ax* + Af + Az\ 
the formula reduces to 

(18) h 2 (x Ax + YAy + ZAz — -^ muA = , 

where X, Y, Z are to be regarded as constant. If $ (Pdp) represents the sum 
of the moments of the impulsive forces, and we regard P as constant, we have 

(19) £ [> (PAp) - 2 ( y mw 2 ) ] = . 

The expressions affected by 8 in these formulae have a greater value than 
they would receive from any other changes of velocity consistent with the 
constraints of the system. 



Deduction of other Properties of Motion. 

The principles which have been established furnish a convenient point of 
departure for the demonstration of various properties of motion relating to 
maxima and minima. We may obtain several such properties by considering 
how the accelerations of a system, at a given instant, will be modified by 
changes of the forces or of the constraints to which the system is subject. 
Let us suppose that the forces X, Y, Z of a system receive the increments 
X, Y, Z', in consequence of which, and of certain additional constraints, 
which do not produce any discontinuity in the velocities, the components of 
acceleration x, y, z receive the increments x', y 1 , z'. The expression 

(20) X [(X+ X')(x + tf)(Y+ Y) (y + y') + (Z+ Z') (z + z>) 

- i- m {(x + if + (y + 'if? + (z + z') 2 }] 

will be the greatest possible for any values of x 1 , y", z' consistent with the con- 
straints. But this expression may be divided into three parts, 
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(21) 2 [(X + X') x + ( Y + T) y+(Z+ Z>) z- \- m (x 2 + f + *)] , 

(22) 2 [Xtf + Zy + Zz' — m (##' + yj/ + zz')] , 
and 

(23) 2 HxV + Y'y' + Z'z' — -|- m (x' 2 +f+ z' 2 )~\ . 

The first part is evidently constant with reference to variations of xf, y", z\ 
and may, therefore, be neglected. With respect to the second part, we observe 
that by the general formula of motion we have 

2 \Xtix + Yhy + Zhz — m (xSx + yhy + z$z)~] = 
for all values of Sx, hy, hz which are possible and reversible before the addi- 
tion of the new constraints. But values proportional to ixf, y", z', and of the 
same sign, are evidently consistent with the original constraints, and when 
the components of acceleration are altered to x + x', y + y', z + z', variations 
of these quantities proportional to and of the same sign as — xf, — /, — z' are 
evidently consistent with the original constraints. Now, if these latter varia- 
tions were not possible before the accelerations were modified by the addition 
of the new forces and constraints, it must be that some constraint was then 
operative which afterwards ceased to be so. The expression (22) will, there- 
fore, be equal to zero, provided only that all the constraints which were ope- 
rative before the addition of the new forces and constraints, remain operative 
afterwards.* With this limitation, therefore, the expression (23) must have 
the greatest value consistent with the constraints. This principle may be ex- 
pressed without reference to rectangular coordinates. If we write u' for the 
relative acceleration due to the additional forces and constraints, we have 

u a — 'xf 1 + f + z' 2 , 
and expression (23) reduces to 

(24) 2 (Xaf + Yjf + Zz' - y miA . 

If the sum of the moments of the additional forces which are considered 

is represented by $(Qdq), (the a representing quantities determined by the 

configuration of the system,) we have 

X(X'x + Y'y+Z>z) = $(Qq). 
d 2 q 

We may distinguish the values of -~ immediately before and immediately 

*As an illustration of the significance of this limitation, we may consider the condition afforded by the 
impenetrability of two bodies in contact. Let us suppose that if subject only to the original forces and con- 
straints they would continue in contact, but that, under the influence of the additional forces and constraints, 
the contact will cease. The impenetrability of the bodies then ceases to be operative as a constraint. Such 
cases form an exception to the principle which is to be established. But there are no exceptions when all the 
original constraints are expressed by equations. 
17 
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after the application of the additional forces and constraints by the expres- 
sions q, and q + q'. With this understanding, we have, by differentiation of 
the preceding equation, 

2 \_X'x + Yy + Z'z + X' (x + '*) + Y (y + y) + Z' (i + *')] 

= » IQg + Q (g. + <?')] ; 

whence it appears that 2 (JV + Yy + Z'z') differs from $» (Qq 1 ) only by quan- 
tities which are independent of the relative accelerations due to the additional 
forces and restraints. It follows that these relative accelerations are such 
as to make 
(25) »{Qfl-2(±-mu?) 

a maximum. 

It will be observed that the condition which determines these relative 
accelerations is of precisely the same form as that which determines absolute 
accelerations. 

An important case is that in which new constraints are added but no new 
forces. The relative accelerations are determined in this case by the condition 

that 2 (— mvr\ is a minimum. In any case of motion, in which finite forces 

do not act at points, lines or surfaces, we may first calculate the accelerations 
which would be produced if there were no constraints, and then determine 

the relative accelerations due to the constraints by the condition that 2 (-»- mvf 2 ) 

is a minimum. This is Gauss's principle of least constraint* 

Again, in any case of motion, we may suppose u to denote the accelera- 
tion which would be produced by the constraints alone, and v! the relative 
acceleration produced by the forces ; we then have 

2 [m {xx 1 + yy' + zz')] = » 
whence, if we write u" for the resultant or actual acceleration, 

2 f -«r- mu 2 J + 2 f-g- mir\ = 2 (y mu"*\ . 

Moreover, differentiating (25), we obtain 

$ {Qh'f) - 2 [m (*&? + y# + z'hi')-] - , 

*This principle may be derived very directly from the general formula (6), or vice versa, for 2 ( -5- mu' 1 ) 
may be put in the form 

the variation of which, with the sign changed, is identical with the first member of (6). 
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whence, since ha", Saf, hy', hz' may have values proportional to q', %', y 1 , z', 

These relations are similar to those which exist with respect to vis viva and 
impulsive forces. 

Particular Equations of Motion. 

From the general formula (12), we may easily obtain particular equations 
which will express the laws of motion in a very general form. 

Let <?«! , da> 2 , etc., be infinitesimals (not necessarily complete differentials) 
the values of which are independent, and by means of which we can perfectly 
define any infinitesimal change in the configuration of the system ; and let 

da>i • dco 2 

o 1= — , ", = -, etc., 
where da 1 , da 2 are to be determined by the change in the configuration in the 
interval of time dt; and let 

d<Oi ■• dco 2 

° 1 = ^' ° 2 =*-' etc - 

Also let TJ = 2 ( y mu j 

It is evident that Z7can be expressed in terms of a lf o 2 , etc., u n o 2 , etc., 
and the quantities which express the configuration of the system, and that 
(since § is used to denote a variation which does not affect the configuration 
or the velocities) , su= ^ Sa dU^ ^ 

Moreover, since the quantities p in the general formula are entirely deter- 
mined by the configuration of the system 

^ = 5; 0i + ^" 2 + etc -' 

where ■— denotes the ratio of simultaneous values of dp and d(d t , when da> 2 
dw t 

etc., are equal to zero, and -j-- , etc., are to be interpreted on the same prin- 
ciple. Multiplying by F, and taking the sum with respect to the several 
forces, we have 

$ (Pf) = Sljih. , + ^2« 2 + etc., 

where ft = f (p £) , a. = »(p£), ete. 

If we differentiate with respect to t, and take the variation denoted by 3, 

we obtain 

£ (P&jp) = Q.M + H 2 fe 2 + etc. 
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The general formula (12) is thus reduced to the form 

< 26 > («'-£) «">+ ("*-£> ^+ ete -=°- 

If the forces have a potential V, we may write 

/-V7\ f dV dTJ\ ,„•• , (dV dU\ .•• , , 

< 27 ) w - «d ^ + Ur - & ^ + etc - 

where -r— denotes the ratio of dV and d^x when io 2 , etc., have the value zero, 

and the analogous expressions are to be interpreted on the same principle. 

If the variations &o x , 5o 2 , etc., are capable both of positive and of nega- 
tive values, we must have 

(28) g=a, §=n 2 , etc., 

(29) d ^ = d ^ d S = d -^ etc. 

^ dtoi do)i ' dcb 2 d(o 2 ' 

To illustrate the use of these equations in a case in which rfox, du 2 , etc., 
are not exact differentials, we may apply them to the problem of the rotation 
of a rigid body of which one point is fixed. If du,_ , du 2 , do 3 denote infinitesi- 
mal rotations about the principal axes which pass through the fixed point, 
fl lt fl 2 , Sli will denote the moments of the impressed forces about these axes, 
and the value of Cwill be given by the formula 

2TJ- (a + b + c)(ol + ol + 4) 2 -(c4 + ig + ol)(aJi + M + co*) 

-\- 2 (b — c) o 2 o B oj -f- 2 (c — a) 030x6)2 + 2 (a — 5) O!0 2 o 3 

-f (b + c) of + (c + a) o* + (a + 5) oi , 
where a, b, and <? are constants, a -\- b, b -{- c, c + a being the moments of 
inertia about the three axes. Hence, 

-pr = (b— c) o 2 o 3 + (J + c) ox, -^ = (<?— • a) o 3 Ox + (c + a) o 2 , 

dot = ( a —b) ^i«2 + (a + J) o,; 

and the equations of motion are 

•• _ (c — 6) 0)2^3 + i?x 

(a—e) wMi 4- ^ 2 

"* = ^ — ' 

(b — a) ojlWz -f- £ 3 
b -f- a 



